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Abstract. A system of simplified equations is proposed to govern the feedback interactions of large-scale flows 
present in laminar-turbulent patterns of transitional wall-bounded flows with small-scale Reynolds stresses gen¬ 
erated by the self-sustainment process of turbulence, itself modeled using an extension of Waleffe’s approach . 
This text is the English adaptation of [18] . It is here complemented by an annex giving the detailed expression of 
the model alluded to above, absent from the French version due to space limitations. 


1 Context 

Flows controlled by viscous effects close to solid walls experience a subcritical transition to turbulence. 
They stay linearly stable up to values of the Reynolds number R high enough to enable nontrivial regimes 
arising from the nonlinearity of Navier-Stokes equations, in competition with the laminar regime. Two 
examples with base flows stable against infinitesimal perturbations for all R are representative of this 
situation: Hagen-Poiseuille flow in a cylindrical tube submitted to a pressure gradient, and plane Couette 
flow (PCF), the simple shear flow developing between two parallel plates in relative translation. Over a 
limited interval in R called the transitional regime, this competition manifests itself in the form of pockets 
filled with small-scale turbulent flow {puffs, spots) co¬ 
existing with laminar flow. Specificities of this situation 
are review from an experimental, numerical, and theoret¬ 
ical point of vue in [T|. The transitional regime of PCF 
and other plane flows is characterized by the presence 
of wide oblique bands, alternatively laminar and turbu¬ 
lent [2], the understanding of which is made difficult by 
the existence of two spatiotemporal scales, one large, that 
of the laminar-turbulent pattern, the other small, inter¬ 
nal to the regions agitated by swift turbulent eddies, as Large-scale back-flow around turbulent 

illustrated in the picture on the right. band fragments in PCF; see [3] for details. 

The origin of turbulence at a local scale is believed to be well understood in terms of dynamics inside 
the Minimal Flow Unit (MFU) [4] where nontrivial solutions “far” from the laminar state are maintained 
by a Self-Sustainment Process (SSP) [5]. Within this framework, exploiting the local coherence of pertur¬ 
bations to the base flow, subsequent theoretical approaches have led to an interpretation of the transition 
using deterministic chaos theory typical of systems with a small number of degrees of freedom [^, thus 
skipping all spatiotemporal aspects most appropriate for systems of experimental interest. Local coher¬ 
ence can however serve as a starting point to a modeling strategy more adapted to the description of 
laminar-turbulent patterns present all along the transitional range [U App. Bj. It indeed supports the 
idea that the dependence of the hydrodynamical fields in the wall-normal direction is more or less frozen 
while the dynamics in the complementary directions (one-dimensional in the pipe case, two-dimensional 
in the other cases) has regained most of its freedom. This is what we shall consider in the following, 
restricting ourselves to the case of PCF though the approach has more general breadth. In particular, we 
shall focus on the large-scale flows HEl that are more easily detected in the intervals between turbu¬ 
lent domains, within the framework -and using the tools- of the general theory of pattern formation in 
far-from-equilibrium systems |5]. 
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2 Modeling 

PCF is the flow that develops between two parallel plates in relative translation in the streamwise direction 
denoted x; y is the wall-normal direction; the spanwise direction z completes the coordinate system. The 
components of the perturbation to the base flow itb = Uy/h are u, v, and w; 2U measures the relative 
speed of the plates and 2h the distance between them. The Reynolds number is defined as R = Uh/v, 
where v is the kinematic viscosity of the fluid. We adopt a unit system based on h for distances and h/U 
for times. In such a system, R is numerically equal to Xjv. In order to simplify the treatment of pressure, 
we use a representation in which v and Q = dzU — dxW are the fundamental variables. The equations 
governing the perturbation to the base flow in this formulation are reviewed and discussed in the book 
by Schmid & Henningson [1^, p. 155ff. 

From previous approaches [1] we shall keep in mind that the wall-normal dependence of the hydro- 
dynamic fields is adequately described at lowest order through a projection onto basis functions keeping 
track of the no-slip condition and the continuity condition: 

dxU + dyV + dzW = Q. (I) 

The experimentally observed coherence of velocity perturbations suggests us to take: 

u = UogQ{y) + Uigi{y) + u', v = Vifi{y) + v', w = W^g^iy)+ Wigi{y) + w', 
with [mn]: 

fi=A(l-y'^f, go = B{l-y^), gi=Cy{l-y^), 

where A, R, and C are normalization constants. (Functions obtained via Empirical Orthogonal Decom¬ 
position of laboratory or computer data might preferably be used.) Amplitudes C/oq, Woq, and Vi are 
functions of space (x, z) and time t. Quantities u',w', v' are residues. Equations governing les amplitudes 
are obtained by projecting the Navier-Stokes equations onto the three basis functions retained. The 
continuity condition immediately yields: 


dxUo + dzWQ = Q and dxUi + dzWi = fiVi, 

where /3 plays the role of a wave vector in the wall-normal direction (/3 = 7r/2 for stress-free boundary 
conditions and P = y/S in the no-slip case with the polynomials given above). 

Setting Zo,i = dzUQ,i — 9xkFo,i, one obtains m- 

(A - dtVi = iy{A^- 2p^A + pi) Vi - (qA - f)(dx(l7oVi) + 9.(WoEi)) 

+ r Idxx (UiUo) + dx. (DiWo + C/oVFi) + dz, (VFobFi)] + Mv ,, (2) 

dtZ^ -\- bdxZi + bdzVi = v {A — po) Zq + sq {pxz {Uq — Wq) + {dzz — dxx){UaWQj\ 

+ Si [dxz {Ul - Wi) + {dzz - dxx)iUiWi)] - So [dziUlVi) - 5,(VFiVi))] (3) 

dtZi + bdxZo = V {A — pi) Z-Y y 2si [dxz{UiUo — WiWo) + {dzz — dxx){UiWo + UoWi)] 

-Mdz{UoVY)-dx{WoVY)\+Mz,. (4) 

where A = dxx + dzz is the Laplacian in the plane of the flow. Coefficients present in these equations 
are integrals of products of the basis functions and their derivatives over y € [—1,1] (see annex), the 
values of which are not interesting in themselves. They can be obtained from |12j where the residues 
were also analyzed in view of going beyond the lowest significant order in a systematic and quantitative 
way. In contrast, here we neglect all terms in M that account for the interaction between the retained 
amplitudes and the residues and between the residues themselves. Accordingly, we get a system with 
a structure identical to that of the three-field model studied in ffUa], except for a velocity-vorticity 
treatment replacing the primitive-variable approach, which yields slightly different coefficients especially 
for the viscous terms (see also the remark at the end of the annex). Numerical simulations of the Navier- 
Stokes equations have shown that this representation of the flow at first significant order contains about 
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90% of the perturbation energy all along the range of Reynolds numbers corresponding to the transitional 
regime [HFig. B6]. Here we shall be satisfied with this approximation though we recognize that, while the 
large-scale back-flows are reasonably well rendered at the qualitative level [111 b], numerical simulations 
of the so-truncated system show that the banded laminar-turbulent coexistence is not appropriately 
reproduced [13] (which, by the way, suggested us to return to direct numerical simulations [3]). 


3 Scale separation 


We are mostly interested in the problem of large scale flows empirically known to exist mostly from 
numerical simulations of Navier-Stokes equations imE] (see figure) but also detected in laboratory 
experiments M- According to the approach developed in [TTJ b] these back-flows stem from source terms 
associated with the Reynolds stress tensor involving small-scale velocity fluctuations; the resulting system 
of equations is however not closed |7l8j . Overcoming this limitation is the purpose of the present work, 
in which we proceed to an explicit separation of slow from fast space-time dependence of the velocity 
field. To make this filtering concrete, in addition to the natural coordinates x, z, t supposed to be rapidly 
varying we introduce slow coordinates X, Z, T. (Coordinate y has been eliminated at the projection step.) 
In this multi-scale approach, noting y the set of variables of interest {Uo,i; Vi; Wop}, we set: 

y :=yiX,Z,T)+y{X,Z,T\x,z,t). (5) 


The subcritical context does not allow us to implement the scale separation using some rigorous system¬ 
atic perturbation expansion as can be performed in the supercritical case to account for the dynamics of 
convection cells or Taylor vortices close to their instability threshold [9] . So we assume a kind of absolute 
scale separation that leads us to consider that, when applied to the fields governed by IM ) the differ¬ 
ential operators involving the natural coordinates are just splitted into operators involving the slow and 
fast coordinates separately. In particular, the continuity condition Q yields: 


dxUo + dzWo = 0 and d^Ui + d^Wi = pVi, 
dxUo + dzWo = 0 and dxUi + dzWi =/3Vi. 


Let us first consider the slowly varying large scales. To obtain their governing equations, it suffices 
to substitute ([^ into (HH and just keep the terms involving the lowest order derivatives in the slow 
coordinates. This yields: 


/S’^OtVi + i^piVi = -f { dxUpVi + dzWQVi ) , ( 6 ) 

BtZq + i'PoZq + bdxZi + bdzV i = so (ydxWiVi — dzUiVi ^ , (7) 

dxZi + vpiZi + bdxZp = si (j^xWpVi — dzUpVi ) , (8) 

in which overlined terms on the right hand side correspond to the slowly variable parts functions of 
the coordinates X, Z, T. Underlined terms are those retained in unb] where the quasi-steady regime 
{dr = 0) was considered, while assuming that the other terms on the right hand side did not contribute for 
statistical symmetry reasons [HI a, §2.4]. This assumption is valid in the absence of “z O — z” symmetry 
breaking but cannot be maintained in full generality since it is not verified in the late stage of spot 
growth [5] and a fortiori when an oblique laminar-turbulent pattern is present [2] . 

The problem is therefore to evaluate Reynolds stresses on the right hand side of (|6j|^ involving the 
average of products of small-scale fast fluctuations. To estimate them, we notice that the small-scale 
flows involved are induced by the SSP that, following Waleffe [S], we approach within the framework of 
the MFU approximation (local view point as opposed to the global perspective related to large scales). 
The eight-variable model that he introduced accounts for a simplified dependence of the flow in the (x, z) 
plane, further Fourier analyzed in series truncated beyond their first-harmonic terms. The modes he chose 
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correspond to a specific phase choice that insures strict spatial resonance between the different velocity 
components involved. Here we keep the MFU assumption considering a periodized domain x 2 x 
and the first harmonic approximation with basic wave vectors a = j and 7 = 27 ^However we 
use a more general decomposition of the hydrodynamic fields permitting phase shifts between Fourier 
modes, which in a certain sense allows us to incorporate the effects of translational invariance in the local 
couplings [5] . Defining a velocity potential (f> and a stream function ip from which the velocity component 
derive as u = — d^ip + dx4>, v = — Ap, w = dxip + so that Q = — Aip, we take the most general 
expressions possible: 


iFo = —Uoz + WqX + Ao cosax + Bq sinacc + Cq cos 72: + Dq sin72; 

+ Eq sin ax cos 72: + Fq cos ax cos 72: + Go cos ax sin jz + Hq sin ax sin 7Z, 
El = —U 12; + VFia: + Hi cosaa; + Bi sinaa: + Gi COS72; + Di sin72 

+ El sin ax cos 72: + Fi cos ax cos 72; + Gi cos ax sin 72: + iJi sin ax sin 7Z, 
<l>i = A!i cos ax + B'l sin ax + C'l cos 7Z + D'l sin 72 

+ E'l sin ax cos 72: + cos ax cos 72; + Gi cos ax sin 72: + sin ax sin 7Z. 


Underlined variables are those introduced under a different name in Waleffe’s original eight-equation 
model (Wa97 in the following), with ‘w’ as a subscript. These are Ui = —1, Dq = —77^/7, Bq = A^ja, 

Fq = —i?w. Hi = —C^ja, El = Dw, C'l = Uw/ 7 , and G'l = E^. The set {C/q; Wq\ Hi; fUi} := G, overlined 
on purpose, obviously correspond to the local values of the components of the large scale flows introduced 
previously. The remaining set, {Hg; Hq; ...; G); Eip} := H, contain all the space-varying contributions that 
are periodic at the period of the MFU. Accounting for the small-scale coherence, these amplitudes are 
supposed to be functions of the fast time t, and the slow coordinates X, Z and T. The system governing 
the dynamics at the MFU scale is obtained in the usual way by inserting the expressions of the velocity 
components that derive from Eq i and Ei as indicated above in ([ 2 jQ, and next separating the different 
trigonometric lines in x and 2 :. A set of 28 equations for 28 unknowns is obtained, that possess all 
the features, lift-up, viscous dissipation, quadratic advection nonlinearities, expected from Navier-Stokes 
equations for wall-bounded shear flows within the MFU framework. Formally, it reads: 


f^A + CA = M{lt)A + AfiA,A), (9) 

^U + C'U = Af\A,A). ( 10 ) 


The format of the RNL proceedings [TH] did not allow us to give the full expression of system ( M 
here to be found in the annex. As to subsystem (10), the original article in French only presented the 
equations for Uq and Ui [equations (25) and (27) in the annex]: 


^Uo + lypoUo = hso[27"(C(Di - CiD[) + k\EiH[ + FiG'i - E[Hi - F[Gi)] (11) 

^Ui + upiUi = \jsi[2^\G'iDo - GoD'i) + k^FIGo + E[Ho - EoH[ - UgG))] (12) 


where kP = a^ + . The first equation has no equivalent in Wa97. The second one closely corresponds 

to the first equation of his system p. 891] for = 1 -|- Ui, in the present notations: 


^C/i + i^piUi = l^si[2^^G[Do - ^FoG[]. 

As an example for subsystem we give here only the equation for Dq [equation (|^ in the annex] 
attached to the component of u in cos( 7 z) independent of x, i.e. the main ingredient of turbulent streaks, 
—Uvj/a in Wa97. It reads: 


^Dq -|- 1/KqDq — bG'i — so7(GoIUo + ^(^{BqFq — AqEq)') 

+ Si [7 (GiIU 1 -I- /3G[U 1 ) -I- ^aj(^{BiFi — HiUi) -|- /1^(H(U( — B[F[)'^ 
- \a^P{AiG'i + H;Gi + BiH[ + B[Hi ))], 


( 13 ) 
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where, in addition to already introduced coefficients, Kq = 7 ^ +j5o- As detailed in the annex, equations 
in are rather complicated but, here, the viscous relaxation term and the lift-up term on the left hand 
side are easily identified {C[ correspond to 14^/7 in Wa97). Of the three uniform amplitudes Wq, Ui, 
and Wi corresponding to the contribution of M ipl) ^ in Q for that equation, only C/i, i.e. — 1, is 
present in Wa97. In our notations the equation for U^j is reduced to: 

^£>0 + vkIDq -bC[ = iayso-Boid) - \a-iSiAiEi - \a^j3siAiG\ + -ff3siC[Ui. 


Let us notice that, in addition to Wa97 that contains 8 equations, system dM possess 
two other closed subsystems with 15 equations. The first one is made of {Wa97;AUo; I/q ; £ 1 }, 
where AUq := {Aq ; Eq ; Bi; Fi ■ H[} and the second one of {Wa97 ; AWq ; Wq ; B[} where AWq := 
{Co ; Go ; i7i; £'i; £{}. These systems are closed in the sense of any trajectory initiated in the corre¬ 
sponding subspaces lives in them indefinitely. Each subsystem AUo and AWo is made of amplitudes 
respectively associated with the generators of infinitesimal translations along x for AUo and z for AWo. 
For example a solution translated by 5x along x can be deduced from any solution of Wa97 provided 
that the perturbation is reduced to Aq = CxBq, Eq = —exFo, Bi = —exAi, Fi = txEi, and E[[ = —CxG'i, 
with Cx = aSx and it is easily checked on (11) reduced to ^Uq + vPqUq = \'ysQK^{EiE[[ -\- FiG{), other 


terms being absent by assumption, that such a solution generates no uniform component t/o. In contrast, 
any arbitrary perturbation in AUo will generically induce some back-flow Uq. The same is true with 
perturbations in AWo and Wq. Back-flows U are therefore naturally part of the dynamics of the complete 
28-amphtude generalized Waleffe system. 

Now, if we assume that this system only describes the small-scale (MFU) flow, and that the corre¬ 
sponding local solution can experience slow modulations governed by (6|7| 8[), we see nonlocal back-flows 
appear via the terms in M (U) A , an example being given by equation (13) with the presence of kUo, 
U 1 and VUi on its right hand side. These terms account for the advection in a wide sense of the small 
scales of the flow, i.e. the amplitudes belonging to A. 


4 Discussion 


In a purely phenomenological approach resting on a Ginzburg-Landau appropriate for a sub-critical 
bifurcation, Hayot and Pomeau m introduced a nonlocal feedback between the large scale recirculation 
and the intensity of local turbulence by an integral condition. System ®il) can be understood as a 
differential counterpart of this condition. 

Now, in a phase-dynamics perspective [9], essential in the context of extended systems of experimental 
interest, the introduction of two groups of additional amplitudes AUg and AWg associated with the 
infinitesimal translation generators in the plane of the flow bring a novel opening apt to raise the rigidity 
inherent in the MFU assumption, since the strictly fixed dimensions of the MFU become in some sense 
softened by the possibility of wavelength and orientation modulations via the analogue of Eckhaus and 
zig-zag instabilities in convection [9], but here in the fully nonlinear regime far from any instability 
threshold. As soon as the trivial contribution Ui producing the S-shaped mean turbulent flow prohle is 
no longer uniform mm, spatiotemporally slowly varying back-flows Uq, Wq and Wi are unavoidably 
generated. This observation leads us to think that an important factor -couplings via has been 

pointed out relative to laminar-turbulent coexistence in transitional wall-bounded flows. The argument 
is not limited to system dWt where terms in Af are neglected but, through the adiabatic elimination 
of residues, would work on any more general effective model able to render the oblique turbulent band 
regime in PCF and other wall-bounded plane shear flows. (It can be noticed that oblique bands were 
observed in numerical simulations of the two-dimensional model extending (Hm beyond first significant 
order [H]). 

In a previous work m, starting from a simplified four-equation Waleffe model, we developed a theory 
of the emergence of periodically modulated turbulent intensity. In a purely phenomenological way, we 
allowed the diffusion of an inhomogenous distribution of the corresponding four amplitudes in an arbitrary 
direction of space and showed that the laminar-turbulent alternation could arise from a Turing mechanism 
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when the diffusivity of the different amplitudes were sufficiently far apart. We also called for the derivation 
of a model from the primitive equations, explaining such a turbulent diffusion. The approach presented 
above is proposed as an element of answer along these lines but we also see that the result is not at all 
in the expected form: The advection of small scales by large scale back-flows PT7| appears to be much 
more natural than the diffusion induced by some anisotropic effective turbulent viscosity that does not 
seem to come out easily from the primitive equations. 

Acknowledgements. The author wants to express his gratitude to G. Kawahara and M. Shimizu (Osaka), 
T. Tsukahara (Tokyo), Y. Duguet (LIMSI), R. Monchaux and M. Couliou (ENSTA-ParisTech) for nu¬ 
merous discussions, in particular within the framework of the JSPS-CNRS program transturb “Localized 
turbulent structures in transitional wall-bounded flows,” as well as to K. Seshasyanan who took part in 
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Annex: Generalized Waleffe model 


For more readability, amplitudes appearing in Wa97 are colored in blue (underlined in the main text); 
the set {AUo;17o} is here painted in red, and the set {AWo;VFo} in green. The guess to be introduced in 
equations is repeated here for convenience with the appropriate colors: 


iFg = —Uqz + WoX + Aq cos ax + Bq sin ax -I- Cg cosqz -|- Dq sinyz 

+ Eq sin ax cos jz + Fq cos ax cos jz + Gg cos ax sin jz + Hq sin ax sin jz, 
El = —Uiz + WiX + Ai cosax -I- Bi sin ax -I- Gi cos 72; -b Di sin 72 

-b El sin ax cos 72 -b Fi cos ax cos 72: -b Gi cos ax sin 'yz + Hi sin ax sin 'yz, 
El = cos ax + B[ sin ax + G( cos 'yz + D[ sin jz 

+ E'l sin ax cos 72: -b F[ cos ax cos 72; -b G'l cos ax sin 72: -b sin ax sin 7Z. 


The velocity and wall-normal vorticity components deriving from these fields read: 

Uo = Uo-d,Eo, Wo=Wo + d,Eo, 

Ui = Ui+ d^Ei -d,Ei, Vi = -AEi, Wi=Wi+ d,Ei + d^Ei, ^ 

Zq = dzUo — dxWi) = —AEq, Zi = dzUi — dxWi = —AEi, 

where ijjo ^^nd Ei denote the periodically varying parts of Eq and Ei defined above. All the coefficients 
introduced in model can be numerically evaluated once the shapes of the Galerkin basis functions 

are specified. Values given below are obtained using the polynomials mentioned in the main text HD- 
Linear viscous coefficients and non-normal linear coupling constants readQ 


Po = i, Pi = T’ = 3 
while constants relative to the nonlinear interactions are: 


Pi = f, 




0.3780, b = 


0.9820, 


0= ^ s 
y 22 

^ VM 

*0 14 


0.8802, 
a 0.8299, 


r = ^ « 0.5590, 


-2.9047, 


Si = ^ 

*1 g 


0.6455, 


r - _ vm 
' “ 4 


Sg = ^ « 0.5590, 


5i = 


-1.6771. 


The model involves a = jix and 7 = 2^ j^z as parameters, ix and (.z being the in-plane dimensions of 

the MFU. Accordingly, in addition to the coefficients given previously, we introduce: 


= a^ -b 7^, 


= a/n, 7 = 7/At, T = 7/c 


— 2 -2 
g = a^ -7A 


9 ' = 2a7, 


which will simplify the expressions of some coefficients of the nonlinear terms. Let 9 be the angle between 
the streamwise direction and the in-plane diagonal of the MFU, then r = tan 9, g = cos 29, and g' = sin 29, 
relate to the shape of the MFU. 

Subset d^ of system dM relative to the periodically varying flow components contains three parts. 
The first one, coming from the identification of amplitudes in Eq, introduces constants 


Ka = «^ + Po, Po, 


= a^ -b 7^ 


■Po, 


S2 — Sg -b / 3 si, S3 — So + 2 ,dsi, 


and reads: 


^ Ag -b i2K^Ag — —abBi — Sg (aUgC/g “b \a'y {CqHq — DqEq)) 

— Si (aBiU 1 -b \a'y {CiHi — DiEi) + {GiF^ + DiG'i)) 

+ S2 {aB[Wi + ia/37 {C[H[ - D[E[) - ^72 {G[Fi + D[Gi)) , ( 1 ) 

^ Bq + vR^Bq = abAi -b Sg (aAgUg -b ^a'y (CqGq — EqEq)') 

+ Si (aAit/ 1 -b \a'y (GiGi — DiFi) — \ {CiE'i -b Di 7 L()) 

- S2 {aA'iWi + ia/37 (G(Gi - D'lF'i) + ^7^ {C'lEi + D'lHi)) , ( 2 ) 


^ In the stress-free case one would get /3 = 7r/2, po = 0, pi = 0^, pi = jd'^, b = lli/2, and b = 7r/2\/2 [111 a]. 
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^ C'o + — —bD[ — So {iDoWo + |a7 {BqGq — AqHq)^ 

- Si (7^?iWi + + ia7 (BiGi - A^H^) + iA[H[ - B[G[) 


+ \a^P (AiF; + A'^Fi + BiE[ + B[Ei)) , (3) 

3^ Do + = bC[ + So (7 GoTFo + {BqFo — AoFg)) 

+ Si (7G1W1 + p-fC[Ui + \a-i (FiFi - AiFi) + \a 0 ^-i {A[E[ - B[F[) 

- \a^P (AiG'i + ^'iGi + BiH[ + B[H^)) , (4) 

^^Eo + vKa-yEo = abFi - ^bH[ + So {aFoUo - iHoWo - a^gAoDo) 

+ Si ^aFiGi - 7F1VP1 - l 3 -fH[Ui - ajgAiDi - \g^ PB'^Ci'^ 

- as 2 F[Wi - \g''^soBiC[ +0/37(0:^50 + gl 3 si)A[D[, (5) 

3^+0 + viil^^Fo = -abEi - 'ybG[ - sq (oFqGo + "fGoWo - ajgBoDo) 

— Si ^oFiGi + 7G1I+1 + PjG'iUi — a^gBiDi + ^g'^f 3 GiA'i'^ 

+ as2E[Wi - \g'^K^S3AiC[ - a/37(o^so + g/ 3 si)B[D[, (6) 

^ Go + vkI^^Go = ^bF[ - abHi + sq {^FqWo - aHoUo - a^gBoCo) 

+ Si (7F1TF1 - aHiUi + /37Fi'l7i - a-fgBiGi - y'^n^PDiA'^ 

Fas 2 H[Wi — y K^sgAiF^ + 0/37(0^50 + /// 3si)i?(^G(, (7) 

^ Flo + ^Q-y^o = otbGi + 7&F( + So (7F0VP0 + otGoUo + cx.'ygAoGo) 

+ Si ^7^1 IFi + oGiG 1 + p^E[Ui + a^gAiCi — y"^I 3 B[D^ 

- 0S2Gyl+i - y'^K^soBiD'^ - a/37(o^so + 5/3si)AyG(. (8) 


A second set stems from amplitudes involved in Fi, introducing constants 

K^=o2 +P1 , k 2 ^ 72 +Pi, 1^2 _^ = q.2 _|_ ^2 _|_ ^ S4 = Si+2/3si. 

It reads: 

^ Ai + vk\Ai = —abBo + si (aj{DoEi + DiEq — Gpifi — GiHo) — 2a{BoU 1 + BiUo) 


- y{GoF[ + DoG\)) + S4 {aB[Wo - y{C[Fo + F^Gq)) , ( 9 ) 

^ i?i + vk^Bi = abAo + si (o7(GoGi + GiGo — DoFi — DiFq) + 2a{AoUi + AiUq) 

- Pi^GoE[ + DoH[)) - S4 (oA'iI+o + yiD[Ho + G[Eo)) , (10) 

Gi + = Si (o7(Aoi7i + AiiJo - FoGi - BiGo) - 27 (FoWi + DiWo) 

- o2/3(AoF{ + BoE[)) - S4 {-fD[Uo + y^iA^Fo + B[Eo)) , (11) 

^ Fi + yOi = Si (o7(SoFi + FiFo - AqFi - AiFq) + 27 (GoWi + Gil+o) 

- ay{AoG[ + BoH[)) + 54 {jG[Uo - y^{A[Go + B[Ho )), (12) 

^ Fi + lynlyi = abFo + si (2d2o(FoGi + FiUo) - 2f-f{HoWi + HiWo) 

- 2o75(AoFi + AiFo)) - y\'^S3 {BoC[ + GoF() 

- o(si + 2a^Psi)F[Wo - l{si + 2fPsMUo, (13) 

FFi+ vnly = -abEo - si {2a^aiEoUi + E 3 U 0 ) + 2727 (Goiri + GiWo) 

+ ‘2oi"fg{BoDi + BiDo)) — y k^S4(AoG( + A^Go) 

+ o(si + 2ays,)E[Wo - 7(si + 2f /3si)G;t7o, (14) 
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^ Gi + = -abHo - Si {2a^a{H,Uo + H^U^) - 2^7 {F^W^ + F^W^) 

+ 2a7g(i3oGi -|- BiCq)) — \g' {A^D'-^ + A'iDq) 

+ 7(^1 + 27^/3si)F{(7o + q;(si -I- 2a^ [3 si)F[[Wq, 

^ iJi + vkI^Hi = abGo + Si {2a^aiGoUi + GiUo) + 2fj{EoWi + EJVo) 

+ 2ajg{AQGi + ^iGq)) — \g' sh{BqD'i F B[Dq) 

— Q!(si -I- 2d;^/3si)Gi I'Fo + 7('5i F 2^'^ j3si)E'iU q. 


The third and last set governs amplitudes involved in <?i and introduces^ 




2 

k'^ = + 20-^ F pi, k!^ = 7 ^- 1 - 2 / 3 ^ 7 ^ +Pi, 

Ca = a{f F 2j3r F q), da = F a^q), 
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Pa-y = + 7 ^ + / 5 ^ 

2^2 


= 7 (r + 2/3r + 7 ^g), dry = \a^T ^(r + 7 ^g) 
Ck = {f F 2/3r -I- K^q). 


= k‘^F2I3^k^ Fpi. 

, o., 2 / 

2/'^, 2 , 


Sa = ^T{K‘^{a'^q -I- r) -I- 2a^(3r), 


e - ir-i 
~ 2' 


(K^( 7 ^g -I- r) -I- 2'^'^(dr), 


It reads: 


(15) 


(16) 


pli-^A'y F vn'^A'y = -t0{G^Fi F GiFq + D^Gy F DyG0 

F daiG[Ho - D[Eo) F ea{DoE[ - GoH[) - CaB'JJo, (17) 

pliB[F vk'^B[ = 7V(Go£;i + GiGq + DoHy F DyHo) 

F da{D'iEQ — G[Go) F CaiGoG'i — DqF[) F CaA'iUo, (18) 

F vk'^C'i = a^r{AQEi F AiFq F BqEi F BiEq) 

F d^ {B[Go - A[Ho) F e0AoH[ - BoG[) - c^D[Wo, (19) 

7^7S F vk'^D'i = a^r(AoGi -I- 2 I 1 G 0 + BqFIi F BiHq) 

F d^ {A[Eo - B[Fo) F e^BoFi - A^E[) + c^G(TTo, (20) 

dl-yFt F uk'^^E[ = y'^n\{BoGy F B,Go) 

F c« {{aFiUo - jH[Wo) F y{j^AoD[ - a^A'M) , (21) 

dUm Fi + F AiGo) 

- c, {{aE[Uo F jG[Wo) F y'ij^BoD[ - a^B[Do)) , (22) 

m'tS G[ F vn'^^G'y = y^K\{AoDy F A,Do) 

F c, {{jFiWo - aH[Uo) F y^BoG', - a^B[Co)) , (23) 

piy H[ F = y^t^riBoDy F ByDo) 

F c, yE[Wo F aG[Uo) F y{a^A[Go - 7 "^oG()) . (24) 

Finally, the equations governing U, subset (10) of system (9,10) in the main text explicitly read: 

0UoF lypoUo = yso[2-f0G[Di - GiD() + K0EyH[ F FiG[ - E[Hy - F(Gi)], (25) 

^ Wo + i^poWo = iaso[2a"(7liS( - Wi^i) + k^{E[Fi F GiH[ - EyF[ - G'yHi)], (26) 

^ G 1 -|- vpiU 1 = i7Si[27^(G(Ilo — Goi7() -I- k^(F’(Go + E[Hq — EqF[[ — FqGi)], (27) 

Wi + 12 P 1 W 1 = yy2a0AoB[ - A'yBo) F k0E[Fo F GoH[ - EoFl - G[Ho)]. (28) 


Checking that energy is conserved by quadratic terms is not as easy as for previous Galerkin mod¬ 
els, due to the uses of the velocity-vorticity formulation that eliminates pressure beforehand. Galerkin 

^ In the stress-free for which pi = one gets k'^ = pi: 7 ~ 7*7’ ^0 = gty, see [5]. 
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projection over the chosen basis indeed does not commute with taking the curl of the primitive equa¬ 
tions that involves taking derivatives with respect to the y coordinate. (With trigonometric lines ap¬ 
propriate to stress-free boundary conditions at the plates, it does.) In the present approach, energy 
conservation is hidden in the complicated expressions of the nonlinear interaction coefficients and coef¬ 
ficients y, in factor of the time derivatives in the equations arising from (Pi. It could have been a little 
more conspicuous by making the following changes ya{A'i,B[), {C[,D[) >->• and 

!->• ya-y{E'i,F{,G[,H[), while complicating further the formal expressions of the coeffi¬ 
cients. As a matter of fact, similar changes were also performed by Waleffe to obtain the final expression 
of Wa87 in 0. 

When restricted to the amplitudes involved in Wa97, our system appears to be more generic than 
Waleffe’s model owing to accidental cancellation implied by trigonometric relations in the stress-free case. 
As a matter of fact, in equation 0 for Ai, a term FqG[ is present that has no equivalent in Wa97: There 
is no term in B^V^ in equation (10.6) for Cw, which invalidates the statement in [^, p. 892, middle of left 
column, about feedbacks among the triad {Fo,C(,Ai} (= {73^, Kv, Cw})- (The detailed consequences of 
this observation have however not been scrutinized, nor the stability of putative fixed points of Wa97.) 

Translational invariance properties can be checked by careful inspection of the equations set. For 
example, it is observed that an infinitesimal x-translation not only generates no as stated in the 
main text, but also that the equations for the amplitudes in AUq are, at first order in = ol5x, just 
proportional to the equations for the corresponding amplitudes in the parent solution of Wa97 with the 
same proportionality coefficient e^, and the appropriate signs. 

Finally, when studying the stability of any nontrivial solution to Wa97, one cans split the tangent space 
into four subspaces. The first one corresponds to amplitude perturbations and involves only perturbations 
within Wa97. The second and third ones are respectively related to x and z phase perturbations, namely 
{AUo ; t/o ;-Di} and {AWq ; VFq ; i3(}. We know from the translational invariance property that they 
possess a neutral mode that is isolated as long as one considers the low-dimensional system obtained by 
sticking to the MFU assumption, but at the origin of continuous phase branches when the assumption 
is relaxed [9]. (The presence of Di and B[ in these groups is easily shown not to break the neutrality of 
the phase perturbations with Uq = 0 and Wq = 0.) The last group, {Wi] Hq;Gi-,Gi-, A[-, E[}, does not 
seem to be related to any specific symmetry. 

Going beyond the remarks in §4 and above is left for future work, in particular the search for explicit 
nontrivial solutions to Wa97 with its no-slip coefficients and their general stability properties. 


